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Abstract. We compare transverse and parallel static susceptibilities of in-plane 
uniaxial anisotropic ferromagnetic Heisenberg films calculated in the framework of 
many-body Green's function theory using single-ion anisotropies with the previously 
investigated case of exchange anisotropies. On the basis of the calculated observables 
(easy and hard axes magnetizations and susceptibilities) no significant differences 
are found, i.e. it is not possible to propose an experiment that might decide which 
kind of anisotropy is acting in an actual ferromagnetic film. 
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1. Introduction 

Jensen et al. pQ have measured parallel and tranverse susceptibilities of a bi-layer Co 
film with an in-plane uniaxial anisotropy, and analysed their results with a many- 
body Green's function theory assuming a spin value of S — 1/2. We have generalized 
their work to multilayers and arbitrary spin in Ref . [2] . In both papers an exchange 
anisotropy was used, because it is easier to treat than the single-ion anisotropy. In 
connection with the reorientation of the magnetization of a ferromagnetic film (with 
an out-of-plane anisotropy) as function of the temperature and film thickness we 
have already discussed similarities and differences between single-ion and exchange 
anisotropies P|. In the latter paper the magnetic dipole-dipole interaction was also 
included. 

In the present paper we calculate within the Green's function formalism anisotropic 
in-plane susceptibilities using the single-ion anisotropy, and compare with the re- 
sults of Ref. fH , where the exchange anisotropy was used. Although we have shown 
in Ref. |3| how the single-ion anisotropy can be treated exactly (for any strength 
of the anisotropy) by introducing higher-order Green's functions, the application to 
multilayers and S > 1 is very cumbersome. This is not the case when using, as we 
did in Refs. [3] and [S] and we do in the present paper, an approximate decoupling 
on the level of the lowest-order Green's functions proposed by Anderson and Callen 
jHj, which however is only a good approximation for small anisotropies, as we showed 
in Ref. [T by comparing with 'exact' Quantum Monte Carlo calculations. In keeping 
with Refs. pQ and [2] we do not include the dipole-dipole interaction, because it is 
nearly isotropic for an in-plane situation. 

The paper is organized as follows. In Section 2 we explain the model and establish 
the Green's function formalism. Section 3 displays the numerical results. In Section 
4 we summarize the results and present our conclusions. 

2. The model and the Green's function formalism 

Although the general formalism is rather similar to our previous work we repeat it 
here to make the paper self-contained. 

The Hamiltonian we use in this paper consists of an isotropic Heisenberg ex- 
change interaction with strength J^i between nearest neighbour lattice sites, a second- 
order in-plane single-ion lattice anisotropy with strength K^^, and an external mag- 
netic field B = (B x ,By,B z ): 

n= - 5E J H(W + CT-EW a 
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- ZilB-St + lB+SZ + B'Si). (1) 
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Here the notation S k = S% ± iS% and B ± = B x ± iB y is introduced, where k and 
/ are lattice site indices and < kl > indicates summation over nearest neighbours 
only. The in-plane lattice directions are the x and z-axes. The field B v will be put 
to zero lateron. 

In order to treat the problem for general spin S, we need the following Green's 
functions 

G8TM = «Sf; (S*r(Srr)) w „ , (2) 

where a = (+, — , z) takes care of all directions in space, rj = ±1 refers to the 
anticommutator or commutator Green's functions, respectively, and n > l,m > 
are positive integers, necessary for dealing with higher spin values S. 
The exact equations of motion are 

<*g%7M = A %7 + ((^ «]-; ^) m {sjT))^ (3) 

with the inhomogeneities 

where (...) = Tr{...e~ /3H ) /Tr(e~ l3 ' H ). The equations are given explicitly by 

ij,r) *Ji'7 

TE4(((^; (Sj) m {S]-) n )) - «SJSf ; 
fc 
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+~ E J»(((srst - s k s+y, (s*r(s-r)) 

-^-GSS" + V^.-™ (5) 

After solving these equations the components of the magnetization can be deter- 
mined from the Green's functions via the spectral theorem. A solution is possible 
by establishing a closed system of equations by decoupling the higher-order Green's 
functions on the right-hand sides. Contrary to Ref . j3] , where we proceed to higher- 
order Green's functions, we stay here at the level of the lowest-order equations. For 
the exchange-interaction terms, we use a generalized Tyablikov- (or RPA-) decou- 
pling 

«S?Sf; (S*r(Srr)) v ~ (S?)G%!f + (S%)G%™. (6) 
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The terms from the single-ion anisotropy have to be decoupled differently, because 
an RPA decoupling leads to unphysical results; e.g. for spin S = 1/2, the terms 
due to the single-ion anisotropy do not vanish in RPA, as they should do, because 
in this case J2i K 2 ^{{Sf) 2 ) is a constant and should not influence the equations of 
motion. In the appendix of Ref. |H| we investigated different decoupling schemes 
proposed in the literature, e.g. those of Lines or that of Anderson and Callen 
jH], which should be reasonable for single-ion anisotropics small compared to the 
exchange interaction. We found the Anderson- Callen decoupling to be most ade- 
quate. It consists in implementing the suggestion of Callen PU| to improve the RPA 
by treating the diagonal terms arising from the single-ion anisotropy as well. This 
leads to 



{((S±S* + S?S?);(S*r(S7r))v 



2<5?>(l " PC* + 1) - (S!S?)})G%™. 



(7) 



This term vanishes for S = 1/2 as it should. 

After a Fourier transform to momentum space, one obtains, for a ferromagnetic 
film with N layers, 3N equations of motion for a 3A-dimensional Green's function 
vector G mn : 



(ul-T)G r 



(8) 



where 1 is the 3Ax3iV unit matrix. The Green's function vectors and inhomogeneity 
vectors each consist of N three-dimensional subvectors which are characterized by 
the layer indices i and j 
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The equations of motion are then expressed in terms of these layer vectors, and 



3x3 submatrices r^- of the 3N x 3N matrix T 
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(10) 

After applying the decoupling procedures © and (J7J), the T matrix reduces to a 
band matrix with zeros in the T^- sub- matrices, when j > i + 1 and j < i — 1. The 
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diagonal sub-matrices r# are of size 3x3 and have the form 

1 H? -H? \ 

r« = o -#/ 



where 



2 * 







^ = Zi+iSfiMq-yJ , 
Zi = B* + J iii+ i(Si +1 ) + J i| i_i(5'f_ 1 ) 

+^2(5f)(l - + 1) - (5?5?>] 

flf = S± + (5f)J«(g-70 + JWi<5S.i) + JM-i<5£i> 



(12) 



For a square lattice, to which we restrict ourselves in the present paper, and a 
lattice constant taken to be unity, 7k = 2(cos k x + cos k y ), and q = 4 is the number 
of intra-layer nearest neighbours. 

The 3x3 off-diagonal sub-matrices Tj,- for j = i ± 1 are of the form 



/ 



-Jij\S i ) 




(13) 
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When treating the monolayer, one can use the spectral theorem for calculating 
the components of the magnetization. This was done in Ref. [H] for the case of spin 
S — 1 and an out-of-plane single-ion anisotropy by using the commutator Green's 
functions. In order to obtain sufficient equations it was necessary, to add equations 
coming from the condition that the commutator Green's functions have to be regular 
at u = 0, which we call the regularity conditions. 

The treatment of multilayers is only practicable when using the eigenvector 
method developed in Ref. j^j. The essential features are as follows. One starts 
with a transformation, which diagonalizes the T-matrix of equation JHJ) 



LTR = ft, 



(14) 



where O is a diagonal matrix with eigenvalues uo T (r = 1, ...,3N). For the problem 
above it turns out that there is one eigenvalue equal to zero for each layer, which 
has to be handled appropriately. The transformation matrix R and its inverse 
R _1 = L are obtained from the right eigenvectors of T as columns and from the 
left eigenvectors as rows, respectively. These matrices are normalized to unity: 
RL=LR=1. 
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Multiplying the equation of motion (jHJ) from the left by L and inserting 1=RL 
one finds 

{ul - Sl)LG™ n = LA™ n . (15) 
Defining Q™ n = LG™ n and A™ n = LA™ n one obtains 

(ul - n)Q™ = A™. (16) 

Q'T' n is a vector of new Green's functions, each component r of which has but a 
single pole 

Amn,T 

grnn,r = _J, 

This is the important point because it allows application of the spectral theorem, e.g. 
[TT] . to each component separately. We obtain for the component r of correlation 
vector C mn = LC mn ( where C mn = ((S z ) m {S-) n S a ) with (a = +, -,z)) 

Amn,T -r i Amn,T 

We emphasize that when (77 = —1), the second term of this equation, which is due to 
the anticommutator Green's function, has to be taken into account. This term occurs 
for uj t = and can be simplified by using the relation between anticommutator and 
commutator 

^4.° = A™>\ + 2C mn >° = L {A^_ x + 2C mn ) , (19) 

where the index r = refers to the eigenvector with uj t = 0. 

The term L A™™_ 1 = vanishes due to the fact that the commutator Green's 
function has to be regular at the origin 

hm = 0, (20) 

which leads to the regularity conditions: 

H X A+1™1 + H X A~^™1 + 2H z A z r ^ 1 = 0. (21) 

For details, see Ref. [3]. 
This is equivalent to 

UA™^ = 0, (22) 
because the left eigenvector of the T-matrix with eigenvector zero has the structure 

L oc{H x ,H x ,2H z ), (23) 

what can be seen analytically. For more details concerning the use of the regularity 
conditions, see Refs. [21 El- 
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We mention an alternative method, published in Ref. ^2], of treating zero 
eigenvalues occurring in the equation of motion matrix, which is based on a singular 
value decomposition of this matrix, and where there is no need for the use of the 
anticommutator Green's function. 

The equations for the correlations are obtained by multiplying equation (JTHJ) 
from the left with R and using equation (|2*2*j) : i.e. 

C = R£LA + R L C, (24) 

where £ is a diagonal matrix with matrix elements S%j = 5ij(e^ Ui — for eigenvalues 
uji 7^ 0, and for eigenvalues u)i = 0. 

This set of equations has to be solved self-consistently together with the regu- 
larity conditions (|2ip. This determines the magnetizations and the moments of the 
magnetizations {(S z ) n ) with n = 2S + 1 for the highest moment, S being the spin 
quantum number. For details see Appendix A of Ref. [3] , where an analogous set of 
similar equations is given more explicitly for the case of the out-of-plane situation. 

The susceptibilities with respect to the easy (xzz) an d hard (Xxx) axes are cal- 
culated as differential quotients 

X ZZ = ((S z (B z ))-(S z (0)))/B z 

X* x = (<5"(S-)> - <5»(0)>)/S», (25) 
where the use of B z w = 0.01/ S turns out to be small enough, see also Ref. 0. 

3. Numerical results 

In this section we show numerical results obtained with the single-ion anisotropy in 
comparison with that from the exchange anisotropy, for which the relevant equations 
were derived in Ref. [2] . As the single- ion anisotropy is not active for S = 1/2 
we will show results for S > 1. In an attempt to obtain universal (independent 
of the spin quantum number S) curves, we have scaled the parameters (J, D, B x w) 
entering the Hamiltonian of Ref. as J/S(S + 1) = J, D/S(S + 1) = D (D 
being the strength of the exchange anisotropy), B x ^ z '/S = B x( - Z \ In the present 
paper we use an additional scaling for the strength of the single-ion anisotropy 
K2/{S — 1/2) = K%. This has been proven to be the proper scaling in Ref.jS], 
because it leads to the correct limit lim T _^o(-^2(^)/-^2(0)) = 1, when calculating 
the temperature- dependent anisotropy by minimizing the free energy with respect 
to the equilibrium orientation angle of the magnetization. 
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3.1 The monolayer with arbitrary spin 

In order to compare results obtained with the single-ion anisotropy and with the 
exchange anisotropy we fit the strength of the single-ion anisotropy to K 2 = 5.625 for 
a square lattice spin S=l monolayer such that the easy axis magnetization (S z )/S 
comes as close as possible to the magnetization obtained previously [2] with the 
exchange anisotropy chosen to be D = 5. The exchange interaction parameter is 
J = 100, and a small magnetic field in x-direction is used, B x = 0.01/5, which 
stabilizes the calculation. The comparison is shown in Fig.l. 

Magnetization: single-ion anisotropy (K 2 =5.625) 
exchange anisotropy (D=5) 



s j h h □ single-ion anisotropy 

53 » H x exchange anisotropy 




20 40 60 80 100 120 140 160 180 200 

temperature 



Figure 1: The magnetization (S z ) /Sofa ferromagnetic spin S = 1 Heisenberg mono- 
layer for a square lattice is shown as function of the temperature. Comparison is 
made between Green's function (RPA) calculations using the exchange anisotropy 
(D = 5, crosses) and the single-ion anisotropy (K 2 = 5.625, open squares) with 
Anderson-Callen decoupling. The corresponding results of mean field (MFT) calcu- 
lations are also displayed. Also shown are the quantities 100 * {S X )/(S + 1) for the 
exchange anisotropy and 100 * (S x ) / S for the single-ion anisotropy; the factor 100 
is introduced to make the curves visible. 

It is surprising that the results for the easy axis magnetization (S z ) are very simi- 
lar over the whole temperature range although the physical origin for the anisotropies 
is very different. An analogous result was observed for the out-of plane situation 
discussed in Ref . [3] . The agreement is not so good for the hard axis magnetization, 
which is a constant for the exchange anisotropy for temperatures below the Curie 
temperature, whereas it rises slightly up to the Curie temperature when using the 
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single-ion anisotropy. In Ref. [2] it was shown analytically that the hard axis magne- 
tization is universal for a scaling {S X )/(S + 1) when using the exchange anisotropy. 
For the single-ion anisotropy a scaling (S x ) / S is found to be more appropriate. 
Comparison is made also with the corresponding mean field (MFT) calculations, 
obtained by putting 7^ = in eqn ([12)1 . showing the well known shift to larger 
Curie temperatures (by a factor of about two for the monolayer) due to the missing 
magnon correlations. 

In Figs. 2 and 3 we show the easy and hard axes magnetizations for a monolayer 
with different spin values 5*. Whereas one observes in Fig. 2 a nearly perfect scaling 
for the RPA calculations with the exchange anisotropy (S = 1/2,1, 3/2, 2, 3, 4, 6, 13/2, 
from Ref. [2] ) and a universal Curie temperature Tc{S) for RPA and MFT, this is 
not the case for the corresponding results with the single-ion anisotropy shown for 
5 = 1, 3/2, 4, 5 in Fig.3, although the violation of scaling is not dramatical. 
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+ 
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RPA+MFT: <S Z >/S, 1 00*<S X >/(S+1 ) 
monolayer: S=1 12 ,1 ,3/2,2,3,4,6,1 3/2 




20 40 60f 80 1 004 1 20 1 40 1 60 1 80 200 

JRPA JMFT temperature 



Figure 2: The magnetizations (S z )/S of spin S = 1/2, 1,3/2,2,3,4,6, 13/2 Heisen- 
berg monolayers for a square lattice are shown as functions of the temperature, from 
Ref. 0. Comparison is made between Green's function (RPA) calculations and re- 
sults of mean field theory (MFT), using the exchange anisotropy strength, D = 5. 
Also shown is the hard axis magnetization, which scales to a universal curve when 
using 100 * (S X )/(S + 1), where the factor 100 is introduced to make the curves 
visible. 
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monolayer with single-ion anisotropy 

S=5, 2,3/2,1 




1 60 1 80 200 



temperature 



Figure 3: The magnetizations (S z )/S of ferromagnetic spin S = 1,2,3/2,5 Heisen- 
berg monolayers for a square lattice are shown as functions of the temperature. 
Comparison is made between Green's function (RPA) calculations using the single- 
ion anisotropy strength of K 2 = 5.625, and the corresponding results of mean field 
theory (MFT). Also shown are the quantities 100 * (S x )/S; the factor 100 is intro- 
duced to make the curves visible. There is only an approximate scaling behaviour. 



Turning to the inverse easy and hard axes susceptibilities Xzz an d Xxx we find 
very similar results for the exchange anisotropy and the single-ion anisotropy. In 
particular in the paramagnetic region (T > Xc U rie) one has a curved behaviour as 
function of the temperature for the susceptibilities in RPA due to the presence of 
spin waves, whereas the corresponding MFT calculations show a Curie- Weiss (linear 
in the temperature) behaviour due to missing magnon excitations. One observes a 
slightly less universal behaviour for the results for the single-ion anisotropy, in Fig. 
4 and 5, when comparing with the results of the exchange anisotropy, see Figs. 2 
and 3 of Ref . [2] . This is connected with the fact that using the exchange anisotropy 
one finds universal values for the Curie temperatures Tq PA (S) and T^ IFT (S), which 
is not strictly the case when using the single-ion anisotropy, see Fig. |3JWe were also 
able to show analytically in Ref. j2] that Xxx * S(S + 1) is universal for T < T c 
when using the exchange anisotropy; this is not the case for the single-ion anisotropy. 
The only difference concerns the curves for the not perfectly scaled RPA results for 
X7z- with the exchange anisotropy the curve with the lowest spin value is left from 
the curves with the higher spin values, whereas the inverse is true for the exchange 
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anisotropy, but this is not a very pronounced effect, and does not lead to a significant 
difference between the results for the various anisotropies. 




Figure 4: 'Universal' inverse easy axis susceptibilities Xzz * S(S + 1) of an in-plane 
anisotropic (due to the single-ion anisotropy) ferromagnetic square lattice Heisen- 
berg monolayer as functions of the temperature for spins S — 5,2, 3/2, 1. Compari- 
son is made between Green's function (RPA) and mean field (MFT) calculations. 
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X2 X *S(S+1 ): single-ion anisotropy 
monolayer: S=5, 2, 3/2, 1 




Figure 5: 'Universal' inverse hard axis susceptibilities Xxx * S(S + -0 °f an in-plane 
anisotropic (due to the single-ion anisotropy) ferromagnetic square lattice Heisen- 
berg monolayer as functions of the temperature for spins S = 5,2, 3/2, 1. Compari- 
son is made between Green's function (RPA) and mean field (MFT) calculations. 



3.2 Multilayers at fixed spin S = 1 

In discussing multilayers with the exchange anisotropy we have considered only the 
case of S = 1/2 in Ref. [2|. In order to compare with results from the single- ion 
anisotropy we have to use a larger spin value because S — 1/2 is not active in this 
case. We restrict ourselves to spin S=l in the following. We have performed also 
calculations with S > 1 which scale with respect to the spin in the same way as in 
the monolayer case. 

In Fig. 6 we compare the Curie temperatures for S = 1 multilayers for exchange 
and single- ion anisotropies in RPA and MFT, using for each layer the same param- 
eters as for the monolayer. Remember that the parameters were fixed such that 
the Curie temperatures for both anisotropies coincide for the monolayer. The Curie 
temperatures for the multilayers N = 2, 19 (for N=19 one is already close to the 
bulk limit) are only slightly lower for the single-ion anisotropy than those for the 
exchange anisotropy. 

In Figs. 7 and 8 we compare easy and hard axes inverse susceptibilities calculated 
with single-ion and exchange anisotropy also for the multilayer case. In order to 
avoid cluttering the figures we restrict ourselves to a multilayer with N=9 layers 
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Curie temperatures : MFT and RPA 
for exchange and single-ion anisotropies 
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Figure 6: Curie temperatures of ferromagnetic spin S=l multilayers are shown as 
function of the film thickness foe RPA and MFT using the exchange (open circles) 
and the single-ion (black square) anisotropies. 

and spin S = 1. For N > 9 the corresponding curves would shift only slightly in 
accordance with the saturation of Tc, see Fig. El with increasing film thickness. We 
display only the RPA results for the multilayer (N=9) and compare with the RPA 
monolayer (N=l) result. Again there is no significant difference in the results for 
both anisotropies. We do not plot the corresponding mean field results which are 
shifted to higher temperatures and show in the paramagnetic region only a linear in 
T Curie- Weiss behaviour, whereas the RPA results have curved shapes due to the 
influence of magnon correlations, which are completely absent in MFT. 
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single-ion+exchange anisotropy: Xi*S(S+1 ) 
RPA: N=1 and 9 layers 




300 



temperature 



Figure 7: The inverse easy axis susceptibilities Xzz °f ferromagnetic films in RPA 
for spin S = 1 for a monolayer (N=l) and a multilayer (N=19) as functions of the 
temperature for single-ion and exchange anisotropies. 

4. Summary and conclusions 

We have applied in this paper a many-body Green's function formalim to calculate 
in-plane anisotropic static susceptibilities of ferromagnetic Heisenberg films using 
the single-ion anisotropy, and compared with previous calculations [2] in which an 
exchange anisotropy was used. Although both kinds of anisotropies are of very 
different physical origin, it is possible, by fitting the strengths of the anisotropies 
properly, to obtain nearly identical values for the easy axis magnetizations over the 
complete temperature range for an S = 1 monolayer. Using the parameters obtained 
in this way also for monolayers with higher spin values and for multilayers, we looked 
for differences in the results of calculations with both kinds of anisotropies. 

By using scaled variables we find a fairly universal (independent of the spin quan- 
tum number S) of easy and hard axes magnetizations and inverse susceptibilities. 
Universality is better established for the exchange anisotropy; e.g. we find a univer- 
sal Curie temperature Tc{S) for RPA and MFT. The scaling is not as perfect for 
the single-ion anisotropy, but there are no dramatic deviations, which might lead to 
a distinction of the influence of both anisotropies. The general statement made in 
Ref. that it is sufficient to do a calculation for a particular S and then to apply 
scaling to obtain the results for other spin values, remains valid to a large extent also 
for the use of the single-ion anisotropy. It remains also true that the measurement 
of the hard axis susceptibility gives in principle the possibility to obtain together 
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single-ion+exchange anisotropy: X X 1 X *S(S+1 ) 
RPA for N=1 and 9 layers spin S=1 

400 
300 
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CO 200 
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100 200 300 
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Figure 8: The inverse hard axis susceptibilities Xxx °f ferromagnetic films in RPA 
for spin S — 1 for a monolayer (N=l) and a multilayer (N=19) as functions of the 
temperature for single-ion and exchange anisotropies. 

with a measurement of the Curie temperature information on the parameters of the 
model, the exchange interaction and the anisotropy strengths. One should, however, 
keep in mind that the quantitative results of the present calculations are due to the 
use of a square lattice. They could change significantly by using different lattice 
types and also by layer-dependent exchange interactions and anisotropies. Such cal- 
culations are possible, because the numerical program is written in such a way that 
layer-dependent coupling constants can be used. 

As a general result we state that our investigations up to now have not lead to 
any significant differences for the calculated observables (easy and hard axes magne- 
tizations and susceptibilities) when using on one hand the single-ion anisotropy and 
on the other hand the exchange anisotropy. Therefore it is not possible on the basis 
of our results to propose an experiment that could decide which kind of anisotropy 
is acting in an actual ferromagnetic film. 
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